Classical mechanism for negative magnetoresistance in two dimensions 
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The classical two-dimensional problem of non-interacting 
electrons scattered by a static impurity potential in the pres- 
ence of magnetic field is investigated both analytically and 
numerically. A strong negative magnetoresistance exists in 
such a system, due to freely circling electrons, which are not 
taken into account by the Boltzmann-Drude approach. A 
parabolic magnetoresistance is found at low fields. 

PACS numbers: 05.60. +w, 73.40.-c, 73.50.Jt 



The negative magnetoresistance, i.e. decrease of resis- 
tance in magnetic field, frequently observed in semicon- 
ductors, as well as in metals, remained a mystery for a 
long time, until Altshuler et al 0] explained this phe- 
nomenon by quantum interference effects (weak local- 
ization). Extensive experimental studies of magnetore- 
sistance, mostly in 2D semiconductor structures, have 
revealed that, apparently, there are two distinct types of 
negative magnetoresistance: (i) a small drop of resistivity 
observed at low fields, such that the classical parameter 
[3 = uj c t is small and (ii) a relatively large (up to 50%) de- 
crease of resistivity at (3 >, 1 or even at j3 » 1. Negative 
magnetoresistance of type (ii) is observed before the on- 
set of Shubnikov-De Haas oscillations and may continue 
as a background trend after the oscillations set in (see, 
for example, Refs. [2-5]. In some instances, the small 
low-field dip of type (i) is superimposed on a smooth 
overall decrease of resistivity ||. Most of the studies 
were devoted to the low-field magnetoresistance of type 
(i), which is very well explained by the weak localiza- 
tion correction for non- interacting electrons pi. As to 
the high-field effect (ii), it is not so well understood, and 
is either attributed to the effect of electron-electron in- 
teraction |^| , which was considered theoretically in Refs. 
H , @] , or left without any explanation. 

We recall that the simple Drude approach yealds zero 
magnetoresistance. The Drude conductivity tensor is 
given by: 



where — ne 2 T jm is the zero-field conductivity, n is the 
electron concentration, e and m are the electron charge 
and effective mass respectively, r is the momentum re- 
laxation time, (3 = lo c t, and lo c = eB/mc is the electron 
cyclotron frequency. For the resistivity tensor it follows: 



Pxx = Po = l/oo, Pxy — ft/ (To — B/nec, and therefore 
the longitudinal resistivity is independent of the mag- 
netic field. This result applies to degenerate electrons 
for which the time r, entering Eq. (1) should be taken at 
the Fermi energy (for nondegenerate electrons one should 
take into account the dependence of the scattering time 
t on the electron energy, which, after averaging Eqs. (1) 
over the Boltzmann energy distribution, results in a pos- 
itive magnetoresistance) . 

In this report, we draw attention to a simple classical 
mechanism for negative magnetoresistance. We consider 
non-interacting 2D electrons with a given energy (the 
Fermi energy) scattered by short-range impurity centers 
in the presence of a magnetic field perpendicular to the 
2D plane, and we show that for any type of scatter- 
ing a strong negative magnetoresistance should exist for 
(3 » 1. We perform computer simulations of the elec- 
tron dynamics in such a system and find an excellent 
agreement between the numerical results and a very sim- 
ple theory which is based on previously known results. 
Moreover we show for the first time that the magnetore- 
sistance is parabolic at low fields. 

The problem was first studied in the pioneering work 
of Baskin et al Q , and more recently by Bobylev et al || , 
who considered specifically the 2D Lorentz model (scat- 
tering by hard disks) and derived the main results rele- 
vant for the further discussion. 

The main idea put forward in Refs. [g, || is that, 
except for the case of small (3, the Boltzmann-Drude ap- 
proach does not work, even as a first approximation, be- 
cause of the existence of "circling" electrons, which never 
collide with the short-range scattering centers, the frac- 
tion of such electrons being |j) 

P = exp(-2irR/e) = exp(-27r//3), (2) 

where R = v/oj c is the cyclotron radius, v is the elec- 
tron (Fermi) velocity, and t = vt is the electron mean 
free path. Contrary to the assumption intrinsic to the 
Boltzmann-Drude approach, an electron which happens 
to make one collisionless cycle will stay on its cyclotron 
orbit forever. The behavior of the rest of electrons (the 
"wandering" electrons, in terms of Ref. |^|), whose frac- 
tion is 1 — P, is controlled by the parameter NR 2 , the 
number of scatterers within the cyclotron orbit, N be- 
ing the impurity concentration. For NR 2 >> 1 they 
behave basically as predicted by the Drude theory, with 
an important modification: after a collision with a given 
scatterer there is a probability P that the electron will 
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recollide with the same scatterer without experiencing 
any other collisions. As a result, for (3 >> 1 the electron 
will recollide with the same impurity center many times, 
and its trajectory will have a form of a rosette, sweeping 
a circular area of radius 2R around the impurity cen- 
ter Since the number of impurities inside this area, 
AttNR 2 , is large, eventually the electron will collide with 
one of them, and thus continue its diffusion in the 2D 
plane. As it follows from the results of Ref. |9| , frequent 
recollisons with the same center lead to the isotropization 
of scattering, so that the effective r in Eq. (1) becomes 
field-dependent. This effect is absent if the scattering is 
isotropic. 

At strong fields, when the parameter NR 2 becomes 
small enough, the rosettes around different scatterers do 
not overlap anymore, the colliding electrons become lo- 
calized and give zero contribution to both a xx and o x 
This means that a percolation transition should occur 
The calculated threshold is {NR 2 ) C = 0.36 §. 

Thus, there are two characteristic values of the mag- 
netic field, B\ defined by ui c = 1/r (/3 = 1), and B2 
defined by oj c = vVN {NR 2 = 1). The ratio B 1 /B 2 = 
(Nd 2 ) 1 / 2 << 1, where d is the scattering cross-section, 
is the small parameter of the theory. 

For the simpler case of isotropic scattering and B « 
B 2 {NR 2 » 1), it follows from the results of Ref. |] 
that the conductivity tensor for wandering electrons is 
given simply by the conventional Drude expressions, Eq. 
(1), with an additional factor (1 — P) in both a xx and 
a xy . The circling electrons behave like free electrons with 
an effective concentration nP, giving a zero contribution 
to a xx , but contributing a term Pgq/P = Pnec/B to 
a xy , and this is the reason why the magnetoresistance is 
negative. This role of circling electrons was overlooked 
in Ref. 9j, but was recognized later jl0| (see also Refs. 

Thus, the conductivity tensor is given by: 
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As a consequence, for the resistivity tensor we obtain 



(36) 



Pxx = P0 



1 - P 

l+P 2 /P 2 



Pxy = PoP 
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1 + P 2 /P 2 ' 



(4) 



Formulas equivalent to Eqs. (3,4) were previously ob- 
tained by Baskin and Entin |l2| for scattering by ran- 
domly positioned antidots. The expression for p xx clearly 
exhibits negative magnetoresistance. Since the terms 
P/P 2 and P 2 1 P 2 are small for any P, Eqs. (4) are very 
similar to 



Pxx = Po(l - P), 



Pxy = PoP = , (5) 

nec 



with an accuracy better than than 2% for p xx , and better 
than 4% for p xy . Note, that at low fields Eqs. (4,5) 
predict an exponentially small magnetoresistance. 

Before further discussion, let us present the results of 
our numerical simulations. In our model a point particle 
(electron) with a given absolute value of its velocity, v, 
is scattered by disks of diameter d randomly positioned 
on a plane inside a square box of edge length L (we take 
L/d = 1000 to be sure that L stays more than an or- 
der of magnitude larger than the electronic mean-free 
path). Periodic boundary conditions are imposed at the 
edges of the square box. Both the hard-disk (Lorentz) 
model, which exhibits anisotropic scattering, and a mod- 
ified model with isotropic scattering are studied. To char- 
acterize the coverage, we introduce a dimensionless con- 
centration of scatterers c = nNd 2 /4, which was changed 
from c = 0.025 to c = 0.2. Studies of the percolation 
phenomena are beyond the scope of the present study. 




FIG. 1. Examples of simulated circling and wandering tra- 
jectories for P — 1 (left) and /3 = 4 (right) for a dimensionless 
impurity concentration c = 7riVd 2 /4 = 0.1. The actual frac- 
tion of circling electrons at yd = 1 is very small. 

In the simulation, we first choose an initial electron 
position at random with an initial velocity along the x- 
direction. In a magnetic field perpendicular to the plane 
the electron trajectory is made of successive circular arcs 
of radius R. For each collision, we determine the inter- 
sections of the trajectory with the disk periphery (the 
impact point), which gives us the impact parameter, b, 
and calculate the scattering angle, 0, accordingly. We fol- 
low the electron velocities v x {t) and v v {t) during a time 
t = 20t sufficient to get reliable results for the integral 
below, and calculate the components of the diffusion ten- 
sor by the standard formula: 



Da = - 

13 2 



< Vi(0)vj{t) > dt. 



(6) 



For each value of field and disk concentration we take the 
average over 10 2 independent disk configurations, and, 
for each configuration, over 10 6 independent trials for 
the initial electron position. Of course, at B = 0, the 
trajectories are straight-line segments and D xy should 
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vanish (this provides a nice test for the numerical pre- 
cision). The conductivity tensor being proportional to 
the diffusion tensor, the components of the resistivity pij 
are calculated as Dij/(D 2 X + D 2 y ), with an appropriate 
normalization. For the Lorentz model numerical calcula- 
tions of this type were previously performed Jll| with an 
emphasis on the percolation phenomenon. 

Fig. 1 shows examples of simulated circling and wan- 
dering trajectories for two values of magnetic field. The 
circling electrons give undamped oscillating contributions 
to the velocity correlation functions in Eq. (6), accord- 
ingly, the integral in Eq. (6), strictly speaking, does not 
converge at t = oo, but is an oscillating function of the 
upper limit. An average over these oscillations is per- 
formed. The same result could be obtained if a small 
damping of these oscillations were introduced (due, for 
example, to weak phonon scattering) . 

The numerical results for p xx as a function of (3 for the 
model with isotropic scattering are presented in Fig. 2 
(top). The resistivity is normalized to the Boltzmann- 
Drude zero-field value, po. The thick line is the theo- 
retical curve predicted by Eq. (4). One can see that the 
theoretical and numerical curves are qualitatively similar 
and the quantitative agreement becomes better as c de- 
creases. In the limit c — > the numerical results converge 
to the theoretical curve, as they should p3). 
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FIG. 2. Numerical results for the resistivity as a function 
of P = lo c t for different impurity concentrations, compared to 
the theoretical curve given by Eq. (4) for the isotropic scat- 
tering model (top) and for the Lorentz model with anisotropic 
scattering (bottom). Circles, squares, diamonds, and trian- 
gles correspond to c = 0.025,0.05,0.1,0.2, respectively. The 
continuous and dashed thick lines are the theoretical curves 
in the isotropic and anisotropic cases, respectively, and they 
are depicted in the inset on a larger scale. Note the surprising 
crossings at /3 = 2. 



Note that for finite c the value of the zero- field resistiv- 
ity is higher than the Boltzmann value po. The relative 
correction for small c is proportional to cln(l/c) and is 
due to recollisions with the same impurity, which are not 
accounted for by the Boltzmann equation [Q . Note also, 
that the numerical results for finite c approach the lim- 
iting theoretical curve from above for [3 < 2, and from 
below for (3 > 2. This may be qualitatively explained as 
follows. On the one hand, at small f} the resistivity for 
finite c is higher than the c — > Boltzmann value due to 
the cln(l/c) correction. We have found analytically that 
in magnetic field this correction decreases quadratically 
in /3, thus giving a parabolic magnetoresistance propor- 
tional to c(3 2 ~ 1/NR 2 . On the other hand, at large (3 
we are on the way to the percolation threshold, where 
Pxx (but not PxyO becomes zero. So, obviously, for large 
(3 and finite c the resistivity should be lower than the 
limiting value given by Eq. (4) Jl3[ . 

Fig. 2 (bottom) displays quite similar results obtained 
for the hard disk Lorentz model (anisotropic scattering). 
The theoretical curve (thick dashed line) was calculated 
using the results of Ref. || for the wandering electrons 
and adding the contribution of circling electrons, as ex- 
plained above. In both cases all the numerical curves for 
different c cross the limiting theoretical curve at the same 
point (3 = 2 (within our numerical precision). We have 
no explanation for this surprising finding so far. 

One of the reasons, why the finite c corrections are of 
interest, is that Eq. (4) predicts exponentially small mag- 
netoresistance for small values of (3. Corrections to this 
formula make the magnetic field dependence parabolic 
and thus define the magnetoresistance at low fields. In 
order to isolate the 1/NR 2 terms in magnetoresistance, 
we look at the difference 

r = ( Pxx(B) ~ PxxjO) } , Pxx{B) ~ Pxx(0) \ m 

1 Pxxifi) >nUm ^ Pxx(0) >th ' 1 ' 

Here the second term in the right hand side is the nor- 
malized magnetoresistance, given by Eq. (4) for isotropic 
scattering, or by a similar formula taking care of the mag- 
netic field dependence of r for the anisotropic case. This 
term represents the limit NR 2 — > 0. The first term is the 
normalized magnetoresistance found numerically. 

This difference, as a function of 1/NR 2 is presented 
in Fig. 3 for both the isotropic and anisotropic scat- 
tering and for different values of c. One can see that 
all the calculated points reasonably fit a universal lin- 
ear dependence, which corresponds to a quadratic de- 
pendence on magnetic field. While such a dependence 
for small (3 could be anticipated, it is surprising that it 
persists for quite large values of (3. Thus at low mag- 
netic field (J3 < 1) the normalized resistivity behaves like 
I — 0.15/NR 2 , the slope being deduced from Fig. 3. For 
f3 > 1 the magnetoresistance is well described by Eqs. 
(4,5) or, for anisotropic scattering, by similar formulas, 
which take care of the magnetic field dependence of r. 
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FIG. 3. The difference, 8, between normalized magne- 
toresistances, calculated in the limit NR 2 — > oo, and ob- 
tained numerically (see Eq. (7)) for the isotropic (empty sym- 
bols) and anisotropic (filled symbols) models, as a function 
of the square of the dimensionless field {B/B 2 ) 2 = 1/NR 2 . 
Circles, squares, diamonds, and triangles correspond to 
c = 0.025,0.05,0.1,0.2, respectively. The data fit a linear 
dependence 8 = -0.15/iVi? 2 . 

The classical approach is justified if the number of Lan- 
dau levels below the Fermi energy is large. However, it is 
irrelevant whether one can describe the individual scat- 
tering events classically (when the electron wave length 
is small compared to the size of the scatterer) or one 
needs a quantum-mechanical description (in the opposite 
case). So long as weak localisation corrections may be 
neglected, the differential scattering cross-section, even 
though calculated quantum-mechanically, may be used 
in the framework of a classical transport theory. We also 
remark that an infinite lifetime of the circling electrons 
is certainly an idealization. In reality, even an electron 
whose orbit initialy avoids the scattering centers, expe- 
riences forces which will gradually change its trajectory. 
However, it is clear that, if the impurity potential de- 
creases fast enough compared to the average distance 
between impurities, the basic features of the model will 
remain valid. The situation is quite different for scatter- 
ing by a long-range random potential, which is typical 
for high-mobility 2D semiconductor structures. Classical 
magnetotransport in this case was thoroughly studied in 
Refs. P. @. 

In summary, for short-range impurity scattering in two 
dimensions a strong negative magnetoresistance exists, 
the main features of which are accounted for by a remark- 
ably simple classical picture. There are two species of 
electrons, wandering electrons which are described by the 
conventional Drude theory, and circling electrons which 
are collisionless and contribute only to the transverse con- 
ductivity. With increasing magnetic field the fraction of 
circling electrons increases at the expense of the wander- 
ing ones, leading to negative magnetoresistance. How- 
ever, as we have shown, at low fields the magnetoresis- 
tance is entirely determined by small corrections to this 
picture, which give a parabolic, rather than an exponen- 
tially small negative magnetoresistance. This mechanism 
may provide a purely classical explanation of, at least, 



some part of the experimental data on negative magne- 
toresistance in two-dimensional structures. 
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